An improved single particle potential for transport model simulations of nuclear 

reactions induced by rare isotope beams 
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Taking into account more accurately the isospin dependence of nucleon-nucleon interactions in 
the in-medium many-body force term of the Gogny effective interaction, new expressions for the 
single nucleon potential and the symmetry energy are derived. Effects of both the spin(isospin) 
and the density dependence of nuclear effective interactions on the symmetry potential and the 
symmetry energy are examined. It is shown that they both play a crucial role in determining the 
symmetry potential and the symmetry energy at supra-saturation densities. The improved single 
nucleon potential will be useful for simulating more accurately nuclear reactions induced by rare 
isotope beams within transport models. 
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The rapid progress in conducting nuclear reaction ex- 
periments using rare isotope beams provides a great op- 
portunity to explore the isospin dependence of strong nu- 
clear interactions The latter determines not only 
the structures of and reactions induced by rare isotopes 
but also the Equation of State (EOS) of neutron-rich nu- 
clear matter relevant for understanding many interesting 
astrophysical phenomena [H4I3 . One of the most im- 
portant inputs for simulating nuclear reactions induced 
by rare isotopes is the single nucleon potential, espe- 
cially its isovector part, i.e., the symmetry potential Q- 
However, our current knowledge about the latter is still 
very poor despite the great efforts made in recent years 
by many people. In fact, a number of theoretical ap- 
proaches, both microscopic and phenomenological in na- 
ture, have been used in studying the single nucleon po- 
tential, such as the relativistic Dirac-Brueckner-Hartree- 
Fock (DBHF) [3 [H], the nonrelativistic Brueckner- 
Hartree-Fock (BHF) approaches [20| . the chiral pertur- 
bation theory [2l[, the nonrelativistic mean- field theory 
using various effective interactions [22J , and the relativis- 
tic impulse approximation coupled with the relativistic 
mean-field models [23T - I25T ]. As one expects, the single nu- 
cleon potentials obtained usually depend rather strongly 
on the details of the model nucleon-nucleon interactions 
used in various many-body approaches, especially at high 
densities and/or momenta [26|-|29(- For instance, within 
the Thomas-Fermi model, it was shown recently that 
the symmetry potential is essentially determined by the 
competition between the isospin singlet (isosinglet T=0) 
and isospin triplet (isotriplet T=l) channels of nucleon- 
nucleon interactions [3(j. Moreover, the resulting sym- 
metry potential is found to be very sensitive to vari- 
ous in-medium effects, such as the in-medium effective 
nuclear many-body forces and the short-range tensor 
force due to the in-medium p meson exchange [3ll l32j . 



It is well known that in nonrelativistic models the in- 
medium many-body force effects can be taken into ac- 
count through a density-dependent term in the two-body 
effective interactions, such as in the Skyrme, the M3Y 
and the Gogny forces [33T - |35"| . In relativistic approaches, 
on the other hand, the dressing of the in-medium spinors 
introduces the density dependence in the two-body inter- 
action leading to similar effects as in the nonrelativistic 
approaches [3|| H3] ■ Interestingly, while all effective in- 
teractions are adjusted to reproduce the saturation prop- 
erties of symmetric nuclear matter, the symmetry poten- 
tials/energies calculated are rather model dependent es- 
pecially at high momenta and/or densities. Therefore, 
it is necessary to further study in more detail the in- 
medium many-body force effects on the symmetry po- 
tentials/energies. 

In this work, taking into account the isospin depen- 
dence of the density-dependent term in the Gogny effec- 
tive interaction we derive new expressions for the single 
nucleon potential and the nuclear symmetry energy. By 
comparing with the old ones previously obtained using 
the original Gogny force [U , we investigate effects of 
the spin(isospin) and density dependence of the many- 
body force term on the symmetry potential and the sym- 
metry energy. The improved single particle potential ob- 
tained with the more complete isospin dependence of the 
effective many-body force term will be useful for more 
accurately simulating nuclear reactions induced by rare 
isotope beams within transport models. 

The central part of the original Gongy effective inter- 
action is 1351 
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where W, B, H, M, and \x are five parameters and P a 
and P T are the spin and isospin exchange operators, re- 
spectively. The last term is an effective two-body force 
deduced from a three-body contact force [3314351] . The 
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xq is the spin(isospin)-dependence parameter controlling 
the relative contributions of the isosinglet and isotriplet 
contributions while the a is the density-dependence pa- 
rameter used to mimic in-medium effects of the many- 
body forces. Since the nn, np, and pp interactions are all 
assumed to depend on the same total density p = p n +p p , 
the proper isospin dependence is neglected in the in- 
medium effects of the many-body forces. However, as 
it was pointed out earlier in Refs. [38l - |4^ there is no a 
priori physical justification based on first-principles for 
any of these interactions to have this kind of density de- 
pendence. In fact, Brucckner, Dabrowski, Haensel et al. 
studied in detail the nuclear symmetry energy within the 
Brueckner theory in the early 60's-70's and found already 
a strong dependence of the G-matrix on the respective 
Fermi momenta of neutrons and protons (k n and k p ) in 
isospin asymmetric nuclear matter (45l l4q | . In another 
word, it is physically more reasonable to assume that the 
interaction between protons depends on the proton den- 
sity, and that between neutrons on the neutron density, 
instead of the total density p [38||. For example, if one 
considers the interaction between two neutrons near the 
surface of a rare neutron-rich isotope or in the neutron- 
skin of a heavy nucleus, it is obviously more appropri- 
ate to assume that the neutron-neutron interaction de- 
pends on the local average neutron density. Indeed, the 
idea of using a separate density-dependence for nn, pp 
and np pairs has already been implemented in various 
models to better understand the structure of nuclei with 
large isospin asymmetries. For instance, local effective in- 
teractions with an appropriate density-dependence sepa- 
rately for nn, pp and np pairs have been proposed a long 
time ago by Sprung and Banerjee (47j . Brueckner and 
Dabrowski |^,[46| and Negele [48|]. Unfortunately, to our 
best knowledge, similar kinds of considerations have not 
been applied yet in simulating heavy-ion collisions involv- 
ing rare isotopes. To help remedy the situation, we adopt 
here the idea from nuclear structure studies. Specifically, 
we replace the density-dependent term in Eq. ([1} with 
the following 

V D = t (l + xoP^WM) + Pr,(rj)] a a(r y ) (2) 

where p T (r) denotes the density of nucleon r (neu- 
tron/proton) at the coordinate r. Using this improved 
density-dependent term, we shall present in the follow- 
ing an analytical expression for the single nucleon po- 
tential that can be used as an input for transport model 
simulations of nuclear reactions. Compared to the MDI 
(Momentum-Dependent Interaction) single nucleon po- 
tential derived in Ref. [22j and used in the IBUU04 trans- 
port model [49l ) . we expect the improved version (IMDI) 
will help more accurately simulate nuclear reactions in- 
duced by rare isotope beams. 

The original MDI single nucleon potential was derived 
from the Hartree-Fock approximation using the original 
Gongy effective interaction [22J. For comparisons, it is 
necessary to first recall the MDI single particle poten- 
tial for a nucleon of momentum p moving in asymmetric 



matter of isospin asymmetry S = (p n — p p )/p and density 
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where r — 1/2 ( — 1/2) for neutrons (protons) and t^t': 
a = 4/3 is the density-dependence parameter; f T {r,p) is 
the phase space distribution function at coordinate r and 
momentum p. The parameters B,C t ,t,C t , t ' and A arc 
obtained by fitting the nuclear matter saturation prop- 
erties [22| . The momentum-dependence of the symmetry 
potential steams from the different interaction strength 
parameters C T>T i and C TjT for a nucleon of isospin r inter- 
acting, respectively, with unlike and like nucleons in the 
background fields. More specifically, C un uk e = —103.4 
MeV while C Hke = -11.7 MeV. The parameters A u (x) 
and Ai (x) are respectively 
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The parameter B and a in the MDI single particle po- 
tential are related to the to and a in the Gogny effective 
interaction via to = f — tt^, and a = a + 1, respec- 
tively. The parameter x is related to the spin(isospin)- 
dependence parameter xq via x = (1 + 2xq)/3. The pa- 
rameter xq determines the ratio of contributions of the 
density-dependent term to the total energy in the isospin 
singlet channel (oc (1 + x^)p a+1 ) and triplet channel 
(oc (1 — x )p a+1 ) (for details, see Ref. (35| ) . For instance, 
xq=1 (xq=-1) means that the density-dependent term 
only contributes to the T=0 (T=l) channel. Thus, by 
varying x from 1 to -1, the MDI interaction covers a large 
range of uncertainties coming from the spin(isospin)- 
dependence of the in-medium many-body forces. In fact, 
the different xq or x parameter used in various Skyrme 
and/or Gogny Hartree-Fock calculations (2(| is respon- 
sible for the rather divergent density dependence of the 
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nuclear symmetry energy |30| . However, we emphasize 
here that the parameter x or x does not affect the EOS 
of symmetric nuclear matter because the x(xq) related 
contributions from T=0 and T=l channels cancel out 
exactly, i.e., oc (1 + x )p a+1 + (1 - x )p a+1 = 2p a+1 . 

The potential energy density correspondin g to the im- 
proved density-dependent term of Eq. © is [44| 
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The corresponding contribution to the single nuclcon po- 
tential obtained from taking the partial derivative of 
the potential energy density with respect to the neu- 
tron/proton density is 



U D (p,r) 
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Replacing properly the parameters to, xq, and a used in 
the Gogny force with the B, x, and a used in the MDI 
interaction Q, the complete expression of an improved 
MDI (IMDI) single particle potential can be written as 
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Some of the parameters have to be re-adjusted to fit the 
saturation properties of symmetric nuclear matter and 
the symmetry energy of E sym (po) = 30 MeV at the nor- 
mal nuclear matter density of po = 0.16//m 3 . More 
specifically, the A' u (x) and A\ (x) are respectively 

2B 
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with B=106.35 MeV. In the case of symmetric nuclear 
matter, as one expects, the IMDI single particle poten- 
tial reduces to the original MDI one, i.e., U(p,0,p,r) = 
U'{p, 0,p,r). The symmetry energy corresponding to the 
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FIG. 1: Momentum dependence of the nuclear isoscarlar po- 
tentials from both the MDI interaction and the IMDI inter- 
action at a density of p=0.32 fm -3 . 

To evaluate quantitatively effects of the isospin and 
density dependence of the in-medium many-body force 
term, we compare in the following the symmetry poten- 
tial and the symmetry energy calculated with the MDI 
and the IMDI single nuclcon potentials. In heavy-ion 
reactions induced by rare isotopes, the global reaction 
dynamics is mostly controlled by the isoscalar nuclear 
potential Uq because of its overwhelming strength com- 
pared to the isovector potential U sym . The latter, never- 
theless, determines all the isospin effects that can be ob- 
served by using some delicate experimental observables 
mostly involving ratios and differences of neutrons and 
protons, see, e.g., ref. 0). The Uq and U sym are related 
to the single nucleon potential by the well-known Lane 
relationship |5(j, namely, U n i p « Uq ± U sym 5. Thus, in 
terms of the neutron (U n ) and proton (U p ) single par- 
ticle potentials, the nucleon isoscalar and isovector po- 
tential can be approximated by Uq = (U n + U p )/2 and 
U sym = (U n — U p )/25, respectively. Both the Uq and 
Usym, especially their momentum dependence, influence 
the density dependence of the nuclear symmetry energy 

via us, mm 
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where t(k) is the nuclcon kinetic energy and kp — 
(37r 2 p/2) 1 ' 3 is the Fermi momentum of nucleons in sym- 
metric nuclear matter. While the symmetry energy is 
calculated exactly using Eqs. (QJ and (fTTj) in this work, 
the above relationship is useful for checking the consis- 
tency and understanding the behaviors of the symmetry 
potential and the symmetry energy. We notice here that 
with the parameter x = 1, by design, the Uq and U sym 
obtained using the IMDI are the same as those obtained 
using the MDI and they both reduce to the predictions 
using the original Gogny force (22[. This can be clearly 
seen from Eq. ((5]) as the fourth term is zero with x = 1 
while other terms stay unchanged as in the original MDI 
potential. 




FIG. 2: Momentum dependence of the nuclear symmetry po- 
tentials from both the MDI interaction and the IMDI inter- 
action at the total nucleon density of p=0.16, 0.32, and 0.48 
fin -3 , respectively. 

As one expects, introducing the isospin dependence in 
the density-dependent term of the effective interactions 
does not affect much the isoscalar potential Uq . As an ex- 
ample, shown in Fig. [1] is the Uq as a function of momen- 
tum at twice the normal nuclear matter density with the 
MDI and IMDI single particle potentials. It is seen that 
there is very little difference between the results obtained 
using the IMDI or MDI. On the other hand, there are sig- 
nificant effects on the symmetry potential U sym and con- 
sequently the symmetry energy E sym (p). Shown in Fig. 
[2] is the momentum dependence of the symmetry poten- 
tials at density p=0.16, 0.32, and 0.48 fm~ 3 , respectively. 
Three typical values of the spin(isospin)-dependence pa- 
rameter x — 1,0, and —1 are used. It is seen from the 
left panel that the symmetry potentials with the MDI or 
IMDI are indeed the same using x=l. For the cases with 
x = and —1, one can see from the middle and right 
panels that the symmetry potentials with the IMDI be- 
gin to deviate significantly from the ones with the MDI 
as the density increases. As a result, one expects that the 
symmetry energy will be significantly different at supra- 
saturation densities with the MDI and IMDI potentials 
using x — and — 1. This expectation is confirmed in Fig. 



[3] where the density dependence of the symmetry energy 
is compared using the MDI and IMDI interactions. Be- 
cause the U sym (k) remains unchanged for x=l (see Fig. 
[5]), it is not surprising that the symmetry energy is the 
same with both the MDI and IMDI interactions using 
x = 1. With x = and — 1, however, the symmetry 
energy with the IMDI becomes significantly stiffer com- 
pared to the MDI case. This effect is consistent with the 
variations of the symmetry potential U sym (k) obtained 
using the MDI and IMDI interactions. 
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FIG. 3: Density dependence of the symmetry energies cal- 
culated using the MDI interaction and the IMDI interaction 
with x=l, and -1. 



From all the expressions for the effective interactions, 
the single nucleon potentials and the symmetry poten- 
tials/energies it is clear that effects of the parameter 
x (or mo) depend on the choice of the parameter a (or 
a) originally introduced to mimic the in-medium effects 
of many-body forces. We thus examine next effects of 
this parameter. Firstly, it is important to note that the 
choice of the parameter a and x is not arbitrary. Besides 
the correlation between them, there are existing experi- 
mental constraints that have to be respected, especially 
the incompressibility Kq = %p^(d 2 E /dp 2 ) Po of symmet- 
ric nuclear matter and the slope of the symmetry en- 
ergy L = 3po dEs j^^ \ po at normal density. While both 
the Kq and L still have some uncertainties mostly be- 
cause their extraction from experimental data is model 
dependent, we use here Kq = 210 ± 20 MeV which is 
consistent within error bars with the ones recently used 
in the literature [5lL |52| . For the slope parameter we 
use L = 88 ± 25 MeV extracted from isospin diffusion 
data within the IBUU04 transport model using the MDI 
with a = 4/3 This range of L is also consistent with 
the ones extracted by using other models (53|. However, 
how the use of the IMDI within transport models may af- 
fect the extraction of the parameter L from experimental 
data remains an interesting question to be investigated. 
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Shown in Fig. [4] are correlations between the Ko and L 
calculated with the MDI and IMDI with x = 1, and -1 
and three values of a, i.e., 0x2,3 of | — and f + 35, 

respectively. As already indicated in Fig. [3j calculations 
with the MDI and IMDI result in symmetry energies with 
significantly different slopes at the saturation density ex- 
cept with x = 1 . With the a between | — ^ and | + ^ , 
both the MDI and IMDI with x = as well as the MDI 
with x = — 1 fall into the area constrained jointly by the 
available constraints on the Kg and L. It is worth empha- 
sizing that the Kq and L only constrains the behaviors 
of the EOS and symmetry energy around the saturation 
point, but not at densities far away from the saturation 
point. 
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FIG. 4: Correlations between the Ko and L with the differ- 
ent spin(isospin)-dependence parameter x and the density- 
dependence parameter a. The filled round symbols denote 
the results using the MDI interaction while the open squares 
are from the IMDI with x=l, 0, and -1, respectively. 




FIG. 5: Momentum dependence of the nuclear symmetry 
potentials from the IMDI interaction with different density- 
dependence at a total density of p=0.16, 0.32, and 0.48 fm -3 , 
respectively. 
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FIG. 6: Symmetry energies calculated by the IMDI inter- 
action with different spin(isospin)-dependence and density- 
dependence. 



Shown in Fig. [5] are the symmetry potentials with the 
three different values of the density-dependence parame- 
ter cti,2,3 and x = 1, and —1, respectively. It is seen that 
the variation of the U sym with a is relatively small except 
at high densities with x = — 1. As the density increases, 
the symmetry potential from the IMDI interaction with 
x=l starts to deviate from the MDI ones. Meanwhile, 
the symmetry potential also shifts downwards with the 
larger a values. However, this is not the case for x=-l. In 
this case, the third term in the IMDI interaction disap- 
pears (see Eq.©) and only the fourth term contributes 
to the symmetry potential. Opposite to the case of x=l, 
the U sym moves upwards and becomes more positive with 
increasing a. The corresponding symmetry energies with 
the different x and a parameters are shown in Fig. [51 It 
is seen that the variation of a can alter appreciably the 
high density behavior of the symmetry energy for any 
given values of x. Unlike the parameter x, however, the 
variation of a has negligible effects around and below the 
saturation density. 

In summary, using different density-dependences for 
the like and unlike nucleon pairs within the Gogny ef- 
fective interaction, we derived new expressions for the 
single nucleon potential and the nuclear symmetry en- 
ergy Effects of both the spin(isospin) and the density 
dependence of the nuclear effective interaction are exam- 
ined. It is found that they play a crucial role in deter- 
mining the symmetry potentials and symmetry energies 
at supra-saturation densities. The improved single nu- 
cleon potential will be used to simulate more accurately 
nuclear reactions induced by rare isotope beams within 
transport models. 
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